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Abstract 

We report a search for CP violation in the decay t~ — >■ 7r~_K's (> Ott") Vt using a dataset of 437 
million r lepton pairs, corresponding to an integrated luminosity of 476 fb~^, collected with the 
BABAR detector at the PEP-II asymmetric energy e^e~ storage rings. The CP-violating decay-rate 
asymmetry is determined to be (—0.36 ± 0.23 ± 0.11)% approximately 2.8 standard deviations from 
the Standard Model prediction of (0.36 ± 0.01)%. 

PACS numbers: 13.35.Dx, ll.30.Er 



1 CP violation has been observed only in the K and 12 

2 B meson systems. However, Bigi and Sanda [l| predict u 

3 that, in the Standard Model (SM), the decay of the t lep- u 

4 ton to final states containing a K^^ meson will also have 15 

5 a non-zero decay-rate asymmetry due to CP violation in le 

6 the kaon sector. The decay-rate asymmetry n 

18 

r (r+ ^ TT+i^O Vr) - r (r- ^ -K-Kl Vr) 

~ r (t+ ^ vr+i^o Vr) + r (r- ^ TT-K", v^) 

21 

7 is predicted to be (0.33 ± 0.01) % for decay times compa- 22 

8 rable to the lifetime r^o of the K'^ meson. In a recent pa- 23 

9 per, Grossman and Nir 2] point out that Sanda and Bigi 24 

10 did not include the interference between the amplitudes 25 

11 of intermediate K'l and which is as important as the 26 



pure -fCg amplitude. Therefore the decay-rate asymme- 
try depends on the reconstruction efficiency as a function 
of the -fCg — !> 7r+7r~ decay time. If the selection is fully 
efficient for decay times that are long compared with the 
ii'g lifetime, then the predicted decay-rate asymmetry is 
almost unchanged relative to the prediction of Bigi and 
Sanda [l|, due to a sign error [2]. 

If the measured decay-rate asymmetry shows a signif- 
icant deviation from the SM value then this could be 
evidence for new physics. No evidence for CP violation 
has been found in related studies by EaBAR and Belle in 
_D+ K'l'K^ decays by the Belle collaboration in a 

study of the angular distribution of the decay products in 
T~ — s> Tr^i^Ts Vr decays or by the CLEO collaboration 
i- 
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This paper presents a measurement of Aq using — > 33 
tt~K'^ (> Ott") i/r and charge conjugate decays. The SM 34 
asymmetry is identical for decays with any number of tt'^ 35 
mesons. If there is an asymmetry due to new-physics 36 
dynamics, then the impact of including modes with one 37 
or more tt" mesons may be different. 33 

The analysis uses data recorded by the BaBAR detec- 39 
tor at the PEP-II asymmetric-energy e~^e~ collider, op- 90 
erated at center-of-mass (CM) energies of 10.58 GeV and 91 
10.54 GeV at the SLAG National Accelerator Laboratory. 92 
The BaBAR detector is described in detail in Ref. frj . In 93 
particular, charged kaons and pions are differentiated by 94 
ionization (dE/dx) measurements in the silicon vertex 95 
detector and the drift chamber in combination with an 96 
internally reflecting Gherenkov detector, with identifica- 
tion efficiency greater than 90% for pions and kaons with 
momenta above 1.5GeV/c in the laboratory frame Q. 
The probability of identifying a pion as a charged kaon 
is less than 2%. An electromagnetic calorimeter made 
of cesium iodide crystals provides energy measurements 
for electrons and photons, and an instrumented fiux re- 
turn detector identifies muons [Oj. For momenta above 
1 GeV/c in the laboratory frame, electrons and muons 
are identified with efficiencies of approximately 92% and 
70%, respectively. Based on an integrated luminosity 
of 476 fb~^, the data sample contains approximately 875 
million r leptons. 

Simulated event samples are used to estimate the pu- 
rity of the data sample. The production of t pairs is 
simulated with the KK2F Monte Garlo (MG) event gen- 
erator [13]. Subsequent decays of the r lepton, contin- 
uum qq events (where q = u,d, s, c), and final-state radia- 
tive effects are modeled with Tauola [U], JETSET [l3|, 
and PHOTOS [l^l , respectively. Passage of the particles 
through the detector is simulated by Geant4 14]. 

The T pair is produced back-to-back in the e+e^ GM 
frame. As a result, the decay products of the two t lep- 
tons can be separated from each other by dividing the 
event into two hemispheres - the "signal" hemisphere 
and the "tag" hemisphere - using the event thrust axis 
[lij . The event thrust axis is calculated using all charged 
particles and all photon candidates in the entire event. 
We select events with one prompt track and a 
7r+ 7r~ candidate reconstructed in the signal hemisphere, 
and exactly one oppositely charged prompt track in the 
tag hemisphere. A prompt track is defined to be a track 
with its point of closest approach to the beam spot be- 
ing less than 1.5 cm in the plane transverse to the e~ 
beam axis and less than 2.5 cm in the direction of the 
beam axis. Furthermore, if a pair of tracks is consistent 
with coming from a /iTg or A decay, or from a 7 conver- 
sion after a mass cut and a displaced vertex cut, neither 
track can be a prompt track. The components of mo- 
mentum transverse to the e~ beam axis for each of these 
two prompt tracks must be greater than 0.1 GeV/c in the 
laboratory frame. The event is rejected if the prompt 



track in the signal hemisphere is identified to be com- 
ing from a charged kaon. A candidate is defined as a 
pair of oppositely charged pion candidates with invariant 
mass between 0.488 and 0.508 GeV/c^; furthermore, the 
distance between the beam spot and the 7r+ 7r~ vertex 
must be at least three times its uncertainty (the tt^ tt" 
will be referred to as the "K^ candidate daughters"). To 
reduce backgrounds from non-r-pair events, we require 
that the momentum of the charged particle in the tag 
hemisphere be less than 4 GeV/ c in the GM frame and be 
identified as an electron (e-tag) or a muon (/i-tag). To re- 
duce backgrounds from Bhabha, fi^fi~ , and qq events, we 
require the magnitude of the event thrust to be between 
0.92 and 0.99. 
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FIG. 1: Invariant-mass distributions for the combined e-tag 
and /i-tag samples. The label in each plot indicates the 
reconstructed decay mode (including the charge conjugate 
mode). Points with error bars represent data whereas the 
histograms represent the simulated sample. The histogram 
labeled as "Signal" includes the r~ — > 7r~_K's (> Otv^)ut-, 
residual r" K-K'^{> 0-K°)ur, and ^ n-K°K°iyr 
modes. All selection criteria (including the likelihood ra- 
tio requirement), except the invariant mass (Mrec) criterion, 
have been applied. The vertical lines and arrows indicate the 
Mroc < 1.8 GeV/ selection criterion. 
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Backgrounds from qq events are further reduced by re-123 
jecting events in which the invariant mass Mrec of thei24 
charged particle (assumed to be a pion), the can-125 
didate, and up to three tt" candidates, all in the signali26 
hemisphere, is greater than 1.8GeV/c^ (see Fig. [IJ. Ifi27 
more than three tt*^ candidates are reconstructed in thei2s 
signal hemisphere, the three with invariant masses clos-129 
est to the tt" mass are included in the calculationoo 
of Mrec and the rest are ignored. The 7r° candidatesm 
are constructed from two clusters of energy deposits ini32 
the electromagnetic calorimeter that have no associatedua 
tracks ("neutral clusters"). The energy of each clusteri34 
is required to be greater than 30MeV in the laboratoryus 
frame, and the invariant mass of the two clusters mustne 
be between 0.115 GeV/c^ and 0.150 GeV/c^ The numbero? 
of events in the r~ — >■ tt^ "Att^Vt mode is small andua 
the corresponding invariant mass plot is not included inug 
Fig. [T] 140 

The imperfect agreement between the Mrcc distribu-141 
tions in the data and MC simulation, seen in Fig. [TJ isi42 
attributed to strange resonances that are not includedwa 
in the simulation. The impact of the modeling of thei44 
T decay modes in the MC simulation on the decay-ratews 
asymmetry is found to be small and is included in thewe 
systematic uncertainties. 147 
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FIG. 2: The likelihood ratio yij) used to distinguish r eventsi69 
from qq events (top plot) and the likelihood ratio y{K^) usedi7o 
to select T decays with a K*^ — >■ 7r"'"7r~ (bottom plot). All se-^^j 
lection cuts, except the plotted likelihood ratio requirement, 
have been applied. Points with error bars represent data while 
histograms correspond to simulated events. The histogram la-^" 
beled as "Signal" includes the t~ — >■ tv' K% (>07r°) Vt, resid-^" 
ual r" -J. K-Kl{> Q-k°)v^, and r" -> -K-K^K^Ur modes.i" 
The vertical lines indicate the selection criteria. i76 

177 

122 A likelihood ratio y(r) is used to distinguish r-pairi7s 



events from qq events, and a second likelihood ratio 
y{K^) is used to reduce the background in the sam- 
ple of TT+TT" candidates. The likelihood ratio 
Ui {xi), where i refers to r or A'^, is defined as yi{xi) = 
CI (xl ) / (£| {x,) + wC\ (x1 ) ) where w is the background-to- 
signal ratio estimated from the MC simulation, £f {C\) 
is the likelihood function for signal (background) events, 
and Xi is the set of variables used for likelihood i. Each 
likelihood function is a product of one-dimensional prob- 
ability distribution functions of the variables Xi obtained 
from the MC simulation. For 2/(t), the variables Xi are 
the visible energy (sum of the energies associated with 
all neutral calorimeter clusters and tracks in the event), 
the number of neutral clusters in the tag hemisphere, the 
number of neutral clusters in the signal hemisphere, the 
magnitude of the thrust, and the component of the total 
momentum of the event transverse to the e~ beam axis 
(calculated from all tracks and neutral clusters in both 
hemispheres). The variables used to construct y{K^) are 
the distance from the beam spot to the decay vertex of 
the -fCg candidate in the plane transverse to the e~ beam 
axis, the invariant mass of the K2; candidate daughters, 
the magnitude of the momentum, and the cosine of 
the polar angle of the K^; candidate. The polar angle is 
the angle between the ifg trajectory and the e~ beam 
axis. The cosine of the polar angle discriminates low- 
angle photon conversions from genuine K'^ candidates. 
All kinematic quantities used in the construction of the 
two likelihood ratios, except for thrust, are determined 
in the laboratory frame. Events are selected if 2/(t) > 0.2 
and y{K^) > 0.4 (see Fig. [5]), in order to minimize the 
contamination from background events while maintain- 
ing a high selection efficiency. 

After all selection criteria are applied, a total of 199064 
(140602) candidates are obtained in the e-tag (/i-tag) 
sample, of which there are 99842 (70369) in the t~ sam- 
ple and 99222 (70233) in the t+ sample. 

The sample contains events from two r decay modes, 
T- K- Kl{>QTT°)vr and ^ -K-R^K^Vr, that also 
have mesons in the final state. The decay — ^ 
TT^ K^K'^Vr satisfies the selection criteria if one of the 
neutral kaons decays into 7r+ 7r~ and the other neutral 
kaon decays into 27r° or appears as a meson. 

The selected candidate sample also contains a small 
background component from r decays not containing a 
in the final state, as well as continuum qq (u, d, s 
and c-quark) events. There is no background from BB 
events. 

The numbers of background events of each type are 
estimated from the MC simulation. The accuracy of 
the background estimation is evaluated by measuring 
the ratios of data to simulated event yields in the re- 
gion ?/(t) < 0.1 and y{K'^) < 0.1. A correction factor is 
then applied to the background yield estimated from the 
Monte Carlo simulation in this region. The correction 
factors are determined to be 0.81±0.03 (0.49±0.03) for 
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179 the qq background and 0.9±0.4 (1.0±0.4) for the non- 

180 K^^ T background in the e-tag (/i-tag) samples, respec- 

181 tively. The total numbers of background events are then 

182 estimated to be 1393±79 (1120±65) for decays and 

183 1401±74 (1055±74) for t+ decays in the e-tag (/x-tag) 

184 samples, where all selection criteria (including the re- 

185 quirements on the two likelihood ratios) are applied. The 

186 uncertainties include the statistical uncertainties from 

187 the sizes of the Monte Carlo samples and the uncertain- 

188 ties of the correction factors. The composition of the 

189 sample is given Table HI 

190 After the subtraction of background composed of qlf 

191 and non-/-Cg t decays, the decay-rate asymmetry is mea-^ 

192 sured to be (—0.32 ± 0.23)% for the e-tag sample and^' 

193 (—0.05 ± 0.27)% for the /x-tag sample, where the errors^' 

194 are statistical. 



TABLE I: Breakdown of the sample after all selection crite- 
ria have been applied. The errors of the decay modes with 
A's are dominated by the uncertainties in the branching frac-^^^ 
tions. The background from other r decays and e^e~ — >■ qTp'' 
background are estimated using the data and MC simulation22s 
samples. 229 



TABLE IL Summary of systematic uncertainties in the decay- 
rate asymmetries. 



Source 


Fractions (%) 




e-tag 


/i-tag 




78.7 ±4.0 


78.4 ± 4.0 




4.2 ±0.3 


4.1 ±0.3 


T- ^ TV- K°K° Ur 


15.7 ±3.7 


15.9 ±3.7 


Other background 


1.40 ±0.06 


1.55 ±0.07 



231 
232 
233 
234 
235 
236 
237 



A control sample of — > h^h^h+{> Ott") Vr (exclud-^"" 
ing ATg TT+TT^ decays) in both data and MC simula-^"^ 
tion, where (ft,+ ) represents a negatively (positively) 
charged hadron, is used to confirm that no significant^"^ 
decay-rate asymmetry is induced by the BaBAR detector^" 
or the selection criteria. The control sample is selected by^"' 
requiring that all charged tracks be prompt tracks, which^"'' 
suppresses AT!? contamination due to its displaced decays"? 
vertex. The asymmetries measured in the simulated and^is 
data control samples agree to within the experimental un-249 
certainties of the measurements, which are 0.12% for the^s" 
e-tag and 0.08% for the /i-tag, and include both statisti-251 
cal and systematic components. These errors are taken^s^ 
as systematic uncertainties on the signal asymmetry (see^s^ 
Table m. 254 

Additional studies show no evidence for any charge-255 
dependent biases in the selection criteria. We find^se 
no decay-rate asymmetry in the MC sample of t~ -^^^ 
TT^ATg (> Ott*') decays (no CP violation is modeled in258 
the simulation) where the error on the decay-rate asym-259 
metrics is 0.14% for the e-tag and 0.17% for the /i-tag26o 
events. We vary the selection criteria around their nom-26i 
inal values, and no significant changes in the asymmetry262 





e-tag 


/i-tag 


Detector and selection bias 


0.12% 


0.08% 


Background subtraction 


0.05% 


0.06% 


K°/K° interaction 


0.01% 


0.01% 


Total 


0.13% 


0.10% 



are observed. The decay-rate asymmetry of the back- 
ground events was studied by examining the events re- 
jected by the likelihood ratio criteria and was found to be 
consistent with zero for both data and MC simulation. 

A recent paper [17] suggests that the decay-rate asym- 
metry will be modified due to the different nuclear- 
interaction cross sections of the K° and K° mesons with 
the material in the detector. This effect is not included 
in the MC simulation. A correction to the asymmetry ac- 
counting for this effect is calculated on an event- by- event 
basis using the momentum and polar angle of the 
candidate together with the nuclear-interaction cross sec- 
tions for neutral kaons, which are related by isospin sym- 
metry to the nucleon cross sections [16j. The kaon- 
nucleus cross sections are determined by using the kaon- 
nucleon cross sections and including a nuclear screening 
factor of A'^ '^^, where A is the atomic weight [l^]. The 
correction, which is subtracted from the measured asym- 
metry, is found to be (0.07 ± 0.01)% for both the e-tag 
and the /x-tag samples. The error includes the statistical 
uncertainty in the MC simulation, the uncertainties in 
the kaon-nucleon cross sections [16|, and an uncertainty 
due to the assumption of isospin invariance. The latter 
effect is taken to be 5% by observing that isospin symme- 
try in pion-nucleon cross sections holds to within a few 
percent. The error on the exponent of the atomic weight 
of the nuclear screening factor is 0.003 17[ and its con- 
tribution to the uncertainty in the asymmetry correction 
is negligible. 

The measured decay-rate asymmetries (after correct- 
ing for the difference in neutral kaon nuclear interac- 
tions) are (—0.39 ± 0.23 ± 0.13)% for the e-tag sample 
and (-0.12 ± 0.27 ± 0.10)% for the /i-tag sample, where 
the first error is statistical and the second is system- 
atic. The systematic uncertainties of the e-tag and /i-tag 
results are almost completely uncorrelated. The small 
correlations in the systematic uncertainties for the two 
samples are ignored when the average is computed. The 
weighted average of the two decay-rate asymmetries is 
(-0.27 ±0.18 ±0.08)%. 

The asymmetry measured at this stage still includes 
other T decays with ATg in the final state. Specifically, the 
decay-rate asymmetry is diluted due to t~ — > Vr 
andr^ — ^ K^K^i^r decays. The measured asymmetry 
A is related to the signal asymmetry Ai and the remain- 
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ing background asymmetries A2 and A3 by: 
/1A1+/2A2 + /3A3 



A = 



/l + /2 + /3 
/1-/2 
/l + /2 + fs 



Ac 



where /i, /2, and /s are, respectively, the fractions of 
T- 7r-is:^(> 0'!T°)vr, ^ K-K°{> 07r")i^^, and 

7r~i^°if "1/7. in the total selected sample, shown in 
Table H Within the SM, Ai = -A2 because the K° in 

tt^K'^ (> Ott") z^t- is produced via a K'^, whereas the 
K° in T- K- Kl{>QTT°)vr is produced via a K°. Fur- 
thermore, A3 = in the SM because the asymmetries due 
to the and will cancel each other. Using the rela- 
tions between Ai, A2, and A3, we can compare our result 
with the theoretical prediction by dividing the measured 
decay-rate asymmetry of ^ = (-0.27 ± 0.18 ± 0.08)% by 
(/i-/2)/(/i+/2+/3) = 0.75±0.04 (the correction is iden- 
tical for the e-tag and /U-tag samples). The uncertainty 
on the correction includes the statistical uncertainty and 
uncertainties in the branching fractions. Finally, the 
decay-rate asymmetry for the t" tt~K'^{> Ott'-') 1/^ 
decay for the combined e-tag and /i-tag sample is calcu- 
lated to be Aq = (-0.36 ± 0.23 ± 0.11)%. 

As pointed out by Grossman and Nir, the predicted 
decay-rate asymmetry is affected by the K'^ — > tt+tt^ 
decay time dependence of the event selection efficiency 
Figure [3] shows the relative selection efficiency, de- 
fined as the selection efficiency normalized to unity in the 
range 0.25 < i/r^o < 1.0. In the < t/rj^a < 1 region, 
the relative efficiency is parametrized with the function^" 
(1 — Ae"'^^*"*"-') , where A, B, and to are constants.^" 
In the 1 < t/T^o < 8 region, the relative efficiency is 
parametrized by a second-order polynomial. Both func- 
tions are constrained to unity at t/rj^o = 1. We use 

^ 318 

this parametrization in Eq. (13) of the Grossman and 
Nir paper Q to obtain a multiplicative correction factor 
of 1.08 ± 0.01 for the decay-rate asymmetry, where the 
error is due to the uncertainty in the relative selection 
efficiency. After applying the correction factor, the SM 
decay-rate asymmetry is predicted to be (0.36 ± 0.01)%. 320 

In conclusion, we have performed a search for CP vi-321 
olation using the r~ — )■ Tr~ K'^ {> Ott'^) decay mode.^^^ 
The decay-rate asymmetry is measured for the first time^^^ 
and is found to be (-0.36 ± 0.23 ± 0.11)%. The measure '" 
mcnt is 2.8 standard deviations from the SM prediction, 
of (0.36 ±0.01)%. 
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FIG. 3: The relative selection efficiency as a function of t/Tp^o 
obtained from the Monte Carlo sample. The top plot shows 
the region < t/T^o < 1 and the bottom plot the region 
1 < t/Tj^o < 8. The solid line is the fit to the points in the 
displayed region. The relative efficiency is normalized to be 
unity for the region 0.25 < t/rj^o < 1.0. 
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